SEMI-CLASSICAL STATES 
FOR NON-SELF- ADJOINT 
SCHRODINGER OPERATORS 

E.B. Davies 
March 1998 



Abstract 

We prove that the spectrum of certain non-self-adjoint Schrodinger op- 
erators is unstable in the semi-classical limit h — > 0. Similar results hold 
for a fixed operator in the high energy limit. The method involves the con- 
struction of approximate semi-classical modes of the operator by the JWKB 
method for energies far from the spectrum. 
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1 Introduction 

It is well known that the complex spectrum of many non-self-adjoint differential 
operators is highly unstable under small perturbations 0, |£| |!|; this has been 
investigated in detail for the Rayleigh equation in hydrodynamics in || Ch 4]. One 
way of exploring this fact is by defining the pseudospectrum of such an operator 
H by 

Spec £ (H) := Spec(H) U {z : \\R(z)\\ > e~ 1 } 

where e > and R(z) is the resolvent of H. It is known that Spec £ (H) contains the 
^-neighbourhood of the spectrum, and that it is contained in the e-neighbourhood 
of the numerical range of H. Theorem 1 states that for Schrodinger operators 
with complex potentials Spec £ (if) expands to fill a region U of the complex plane 



much larger than the spectrum in the semi-classical limit h — > 0. More precisely 
we obtain an explicit lower bound on ||i?(z)|| which increases rapidly as h — > for 
all z in the region U defined by 

U ■= {z = 7] 2 + V{a) : 77 G R\{0} and lm(V'(a)) ^ 0}. 

For h = 1 we deduce that for large enough z within a suitable region the resolvent 
norms of Schrodinger operators with complex potentials become very large even 
though z may be far from the spectrum. This had apparently not been noticed in 
other spectral investigations of these operators [p], £|, || [7[, apart from ||, whose 
results are greatly extended and improved in Theorem 2 below. 

Our results have a positive aspect. Our proofs use a JWKB analysis to construct a 
continuous family of approximate eigenstates, which we call semi-classical modes, 
for the operators in question. These modes have complex energies far from the 
spectrum, but could be used to investigate the time evolution of fairly general 
initial states by expanding these in terms of the modes. 



2 The estimates 

For reasons which will become clear in the next section, we consider operators 
somewhat more general than those described in the last section. We assume that 
H acts in L 2 (R) according to the formula 

A 2 -P 

Hf(x) := -h 2 ^-L + V h (x)f(x) 

Where Vh are smooth potentials for all small enough h > which depend contin- 
uously on h, together with their derivatives of all orders. In the applications in 
the next section Vh has an expansion involving fractional powers of h, but this is 
invisible since we treat Vh as it stands, and only need an asymptotic expansion 
involving integer powers of h; this simplification is essential to our solution of the 
problem. Technically we assume that H is some closed extension of the operator 
initially defined on C£°(R). We assume throughout the section that 

z := rf + V h (a) 

where r\ G R\{0} and Im(Vo(a)) 7^ 0. The /i-dependence of z as defined above 
may be eliminated by a uniformity argument spelled out in the next section, but 
we do not focus on this issue here. Our goal is to prove the upper bound 

||^/-2f/||/||/||=0(/l B ) 
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as h — > 0, for all n > 0, where / G C£°(R) depends upon h and n. This immediately 
implies that \\(H — diverges as h —* faster than any negative power of 

h. Although the above equation may be interpreted as stating that such z are 
approximate eigenvalues for small enough h > 0, it does not follow that they are 
close to true eigenvalues, and indeed the examples studied in [Q, [| |, |9|] show that 
there is a strong distinction between the spectrum and pseudospectrum. 

Theorem 1 There exists 5 > and for each n > a positive constant c n and an 
h- dependent function f G C£°(R) such that if < h < 5 then 

\\Hf-zf\\/\\f\\<c n h\ 



Proof The proof is a direct construction. We put f(a+s) := £(s)/(s) for all s G R 
where £ G C^°(R) satisfies £(s) = 1 if \s\ < 5/2 and £(s) = if |s| > 5, and 5 > is 
determined below. We take / to be the smooth but not square integrable JWKB 
function / := exp(— ip) where 



and ip m are defined below. A direct computation shows that 



(2n+2 \ 
E h m K f 
m=0 ) 



where m are given by the formulae 



0i 

02 
03 

?2n+2 



-(V>-i) 2 + ^-* 
V>-i-2V-i$> 

-«) 2 - 



By setting m = for < m < n + 1, we obtain a series of equations which enable 
us to determine all ip m , provided 5 > is small enough. 

The key equation is the complex eikonal equation 

^ l {sf = V h {a + s)-V h {a)-r l 2 
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whose solution is 

V-i(s) := f a (V h (a + t)-V h (a)-,f) in dt 

iVl(a)s 2 ^. o 
477 

We have assumed that Im(V^(a)) 7^ and this implies that Im(V^(a)) 7^ for all 
small enough /i > 0. We choose r\ to be of the same sign as Im(V h '(a)) so that for 
a suitable constant 7 > we have 

7s 2 < Re(^_i(s)) < 37s 2 

for all small enough s and h. We also assume that s and /i are small enough for 
p := (2ip'_ 1 )^ 1 to satisfy a bound of the form \p(s)\ < /3. 

We now force m = for < m < n + 1 by putting 

= pW-W) 

etc. 

We determine the functions uniquely by also imposing ip m (0) = for all < m < n. 
Each of the functions is bounded provided s and h are small enough, and the same 
is true of the remaining functions <f> m . Specifically we assume that for some 5 > 
and constants c m , c' m we have 

\ipm(s)\ < c m 

for < m < n, and 

\<t>m(s)\ <c' m 

for n + 2 < m < In + 2, provided \s\ < 5 and < h < 5 2 . 

In the following calculations dj denote various positive constants, independent of 
h and s. We have 

> 

2 — 



„ l/(«)| 2 d S 

-5/2 

> f S/2 e-^-'-^ds 

J -511 
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-6h- 1 / 2 /2 

r 1 / 2 9 
> / e-^-^h^dt 

7-1/2 

= a 2 /^ 2 . 

We also have 

\\Hf-zfh = \\-h 2 fe-2h 2 f'? + Z{Hf-zf)h 

2n+2 

< h 2 \\ffh+2h 2 \\fe\\ 2 + e m^/ih 

m=n+2 

and need to estimate each of the norms. Since £' has support m {s : 8 /2 < \s\ < 8} , 
we have 

U"f\\l <aJ e-^-^da 

J5/2<|s|<5 

< a 5 e-^ 2 /^. 

In other words H^'/lb decreases at an exponential rate as h — > 0. A similar 
argument applies to ||£'/'||2- Since m is bounded on {s : |s| < 5}, uniformly for 

\h\ < 8 2 , we see that 

UKfWl < a 6 [ 5 \f(s)\ 2 ds 



1-8 

< / Q-^'+^dS 

-8 



< ash 1 / 2 

by an argument similar to that used for / above. Putting the various inequalities 
together we obtain the statement of the theorem. 



3 High Energy Spectrum 

By a change of scale our theorem can be applied to prove the instability of the high 
energy spectrum of a non-self-adjoint Schrodinger operator with complex potential. 
The results in this section extend those of H both by providing greater insight 
into the machanism involved and by obtaining much stronger estimates. Adopting 
quantum mechanical notation we assume that h = 1 and that the operator H 
acting in L 2 (R) is given by 

H := P 2 + J2 c m Q m 

m=l 
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where n is even and the constant c n has positive real and imaginary parts. 
Theorem 2 If z E C satisfies < arg(z) < arg(c n ) and a > then 

\\(H-azin\ 

diverges to infinity faster than any power of a as a — > +oo. 
Proof If u > then the operator H is unitarily equivalent to the operator 

n 

Hl := u- 2 P 2 + £ c mM m Q m 

m=l 

Putting m := cr 1 /™ we obtain 

\\{H-azI)- 1 \\=a- l \\{H 2 -zI)- 1 \\ 

where 

Putting h := ■u~(™+ 2 )/ 2 we have 

n 

H 2 := /i 2 P 2 + c m h 2{n - m)/{n+2) Q m = h 2 P 2 + V h (Q). 

m=l 

This is precisely the form of operator to which Theorem 1 applies. We have Vq(o) = 
c n a n so Im(c n ) > implies that the conditions of Theorem 1 are satisfied for any 
z in the sector 

U := {z : < arg(^) < arg(c n )}. 

This completes the proof, except for a technical point which we now address. In 
Theorem 1 we assumed that z = rf + Vh(a) where Im(V^(a)) 7^ 0, so z is apparently 
dependent on h, with the limit zo := rj 2 + Vo(a) as h — > 0. We rectify this problem 
by fixing z G U and making a and rj depend upon h in such a way that 

z = Vl + V h {a h ) 

where r)h — > 77 and — > a as /i — > 0. We now have to check that all the estimates 
of Section 2 are locally uniform with respect to 77 and a, so that the result we claim 
does indeed follow. 

The method of this paper can be extended to treat certain rotationally invariant 
problems in higher space dimensions. The condition —2 < p(l) in the next theorem 
is included because it is relevant to the existence of a closed extension of the 
operator, by virtue of an application of the theory of sectorial forms. 



n 

= u - 2 -"P 2 + ]T c m u m ~ n Q 

m=l 
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Theorem 3 Let the operator H acting in L 2 (R N ) be some closed extension of the 
operator given by 

n 

Hf(x) := -A/(:r) + ]T c rn \x\^f(x) 

m=l 

for all f in the initial domain C^ >o (R Ar \{0}) ; where c n has positive real and imag- 
inary parts, p{n) > and 

-2 < p(l) < p(2) < . . . < p(n). 
If z e C satisfies < arg(^) < arg(c n ) and a > t/ien 

IK^-^/)- 1 !! 

diverges to infinity faster than any power of a as a — > +oo. 

Proof The difference from Theorem 2 is that after restricting to the usual angular 
momentum sectors the operators act in L 2 (0, oo) and include angular momentum 
terms in the potential. However, it may be seen that the analysis of Theorem 2 
can be extended to operators of the form 

H:=P 2 +jr c m Q p ^ 

m=l 

so the incorporation of the angular momentum terms causes no difficulties. 

Note Since the supports of the test functions used in the proof of the theorem 
are compact and move to infinity, Theorem 3 remains valid if we add a non- 
central potential to H, provided that potential decreases at infinity faster than 
any negative power of \x\. Weaker versions of Theorems 2 and 3 hold if one adjoins 
a potential which decreases more slowly at infinity. 
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